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Abstract 

The /3-decay process of the 6 He halo nucleus into the a + d continuum is studied in a three-body 
model. The 6 He nucleus is described as an a + n + n system in hyperspherical coordinates on a 
Lagrange mesh. The convergence of the Gamow-Teller matrix element requires the knowledge of 
wave functions up to about 30 fm and of hypermomentum components up to I\ = 24. The shape 
and absolute values of the transition probability per time and energy units of a recent experiment 
can be reproduced very well with an appropriate a + d potential. A total transition probability of 
1.6 x 10 -6 s^ 1 is obtained in agreement with that experiment. Halo effects are shown to be very 
important because of a strong cancellation between the internal and halo components of the matrix 
element, as observed in previous studies. The forbidden bound state in the a + d potential is found 
essential to reproduce the order of magnitude of the data. Comments are made on i?-matrix fits. 
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I. INTRODUCTION 


The discovery of light halo nuclei with very large matter radii near the neutron drip 
line |H inspired detailed studies of the structure of these quantum systems. The larce radii 
are interpreted as arising from an extended spatial density of a few neutrons 0, 0. The 
static properties of the halo nuclei do not provide a complete picture of their structure and 
especially of the halo extension. Few observables directly probe the probability density at 
very large distances. 

The f3 decay with emission of a deuteron, also known as f3 delayed deuteron decay, is 
energetically possible for nuclei with a two-neutron separation energy S- 2 n smaller than 3 
MeV. This property is typical of halo nuclei. The measurement of the spectrum shape for 
this decay process offers a unique opportunity of probing halo properties at large distances. 
The (3 decay of 6 He into a and a deuteron has been observed in several experiments B00- 


The branching ratio is much smaller than expected from simple i?-matrix [4], two-body 7], 
and three-body [[§jj models. Various experimental values of this branching ratio have been 
obtained, i.e., (2.8 ± 0.5) x 1(T 6 0, (7.6 ± 0.6) x 10~ 6 0, and (1.9 ± 0.8) x 1(T 6 0, for a 
deuteron cutoff energy of about 350 keV. The latest result |6] is a factor of 4 smaller than 
the result of Ref. 0 which served as reference for most theoretical papers. 

A semi-microscopic study |9] of the process has been able to explain that the low value 
of the branching ratio is the result of a cancellation between the ’’internal” and ’’external” 
parts of the Garnow-Teller matrix element. The overlaps of the 6 He ground state and a + d 
scattering wave functions in the internal (R < 5 fm) and external (R > 5 fm) regions 
have very close magnitudes but opposite signs. It is clear that the external part of the 
Garnow-Teller matrix element reflects properties of the halo structure of the 6 He nucleus. 
An improved microscopic wave function of 6 He confirmed this interpretation [lfij. It was also 
confirmed by a fit in the /?-matrix framework 11] which yields a satisfactory description of 
the deuteron spectrum shape and branching ratio of Ref. 0. A fully microscopic description 
of the (3 decay of the 6 He nucleus to the 6 Li ground state and to the a + d continuum jl2| was 
performed in a dynamical microscopic cluster model with consistent fully antisymmetrized 
wave functions for the initial bound state and the final scattering state. This model provided 
a reasonable agreement with the data of Ref. 0. Without any fitted parameter, those data 
were underestimated by about a factor of 2. Hence, the same microscopic results now 
overestimate the recent data of Ref. 0 by a similar factor. 

Since new data [6] with much better statistics which provide an even lower branching 
ratio are now available, it is timely to reexamine the interpretation of the (3 delayed deuteron 
decay. Improving significantly the microscopic model of Ref. [121] is not yet possible. We 
prefer thus to base our discussion on an a+n+n three-body model. Accurate wave functions 
of 6 He are available in hyperspherical coordinates 0 - A previous calculation based on the 
same model 
data of Ref. 


contains several limitations which led to a significant overestimation of the 
: the calculations were restricted to small values of the hypermomentum, 
K = 0 and 2, and the halo description may not have been sufficiently extended. 

The aim of the present work is the determination of the deuteron spectrum shape and 
branching ratio for the (3 decay of the 6 He halo nucleus into a + d with an accurate treatment 
of the 6 He wave function in an a + n + n three-body cluster model. For the description of the 
structure of the 6 He nucleus, we use the hyperspherical harmonics method on a Lagrange 
which yields an accurate solution in this model. The a + d scattering wave 


mesh jl 

function is treated as factorized into a deuteron wave function and a nucleus-nucleus relative 
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wave function. We will choose several versions of the central interaction potential between 
a and d: a deep Gaussian potential H which fits both the s-wave phase shift and the 
binding energy of the 6 Li ground sta te ( 1.473 MeV), and potentials obtained by folding the 
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a + N potential of Voronchev et al. jl 6 j. For the sake of comparison we will also perform a 


calculation with a repulsive a + d potential which was used in Ref. 

In Sec. II, the formalism of the /3-decay process of the 6 He nucleus into the a+d continuum 
is presented. The potentials and the corresponding three-body hyperspherical and two-body 
scattering wave functions are also described. In Sec. Ill, we discuss the obtained numerical 
results in comparison with experimental data. Finally conclusions are given in Sec. IV. 


II. MODEL 

A. 6 He wave function 


The initial three-body wave function is expressed in hyperspherical coordinates. Particles 
1 and 2 are the nucleons and 3 is the a cluster. A set of Jacobi coordinates for the three 
particles with mass numbers Ai — 1 , A 2 — 1, and A 3 = 4 is defined as 


X = v / AL 2 r 5 y = \fmnyiR; (1) 

where r = r 2 — rq and R = r 3 — l(r 1 + r 2 ). The (dimensionless) reduced masses are given 
by /U 12 = 1/2 and /i(i 2)3 = 4/3. Equations © define six coordinates which are transformed 
to the hyperspherical coordinates by 

p 2 = x 2 + y 2 , a = arctan (y/x), ( 2 ) 


where a varies between 0 and n/2. With the angular variables fl x = (6 x ,ip x ) and £} y = 
(O y ,ip y ), Eqs. © define a set of hyperspherical coordinates. These coordinates are known 
to be well adapted to the three-body Schrodinger equation. 

With the notation G 5 = ( a , fl x , f^), the wave function reads 3 


IJyLSK 

where l x and l y are the orbital momenta associated with the Jacobi coordinates x and y , 
respectively, y/J^sK are hyperspherical harmonics, and Xij y LSK are hyperradial functions. 

The 6 He ground state wave function contains components with total intrinsic spin S = 0 
and 1. The total orbital momentum L is equal to S. Because of the positive parity, l x + l y 
is even and the sums in Eq. © and in the following run over even K values only. 


B. a + d wave function 

For the scattering state, we assume an expression factorized into the deuteron ground- 
state wave function and an a + d scattering wave function derived from a potential model. 
The deuteron spin 1 and positive parity allow S and D components. Here, we neglect the 
small D component of the deuteron. 
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Below, we only need the 1 + component of the a + d scattering function which reads 


+ (E , r, R ) = R)- (4) 

where r and R are here the deuteron and a + d relative coordinates, respectively. The 
spatial part of the deuteron wave function is written as 

M r ) = r -1 u d (r)y 0 o(^x)- (5) 

The spatial part of the a + d s- wave function is factorized as 

ipad(E,R) = R~ 1 u E (R)Yio(n y ). ( 6 ) 

The radial scattering wave function u E has the asymptotic behavior, 

u e (R) -»• cos6 0 (E)F 0 (kR) + sm6 0 (E)G 0 (kR), (7) 

R—> oo 


where k is the wave number of the relative motion, Fq and Go are Coulomb functions, and 
5 0 (E) is the s-wave phase shift at energy E. 

The radial functions u E are calculated with effective a + d potentials. Some among the 
potentials V a d(R) we are using are obtained by folding an a + N potential V aE r(r). They 
are given by the equation 

V ad (R) = (Mr) I V an (\R+ fr|) + V ap (\R-±r\) \ ^ d (r)), (8) 

where the integration is performed over the radial and angular parts of variable r. 


C. Transition probability per time and energy units 


For the /3-decay process 

6 He —> ol T d T e v , 

the transition probability per time and energy units is given by na 


dW 


771 e C 


dE 7 T A vh 2 


G 2 0 f(Q - E)Bgt(E), 


( 9 ) 


( 10 ) 


where m e is the electron mass, v and E are the relative velocity and energy in the center 
of mass system of a and deuteron, and Gp = 2.996 x 10~ 12 is the dimensionless /3-decay 
constant [18(- The Fermi integral f(Q — E ) depends on the kinetic energy Q — E, available 
for the electron and antineutrino. The mass difference Q between initial and final particles 
is 2.03 MeV. 

Between an initial state with isospin T = 1 and a final state with isospin T = 0, Garnow- 
Teller transitions are allowed. The reduced transition probability reads 


BcME) = 12A 2 ^|(»" + (£)|y/t,_^|v|)»+)|^ (11) 

M j =1 

where A = 1.268 ± 0.002 is the ratio of the axial-vector to vector coupling constants 0. 
T^ /+ (i3) is the wave function (J3J) of the final a+d system, and is the wave function @ 
of 6 He. The operators Sj and t 3 are the spin and isospin operators of particle j, respectively. 
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Since the total orbital momentum and parity are conserved, only the l = 0 partial scat¬ 
tering wave contributes. Hence, only the initial L = S = 0 component of 6 He can decay 
to a + d. It is convenient to express the Gamow-Tcller matrix element with the help of an 
effective wave function ||, 


f/’efiK-R) 


= / ’Pd(r)tpS He (r,R)dr. 


( 12 ) 


In this expression, t/’e' H (t*, is the spatial part of the S = 0 component of the 6 He wave 
function. The reduced transition probability can be written as 


■®gt(^') — 6 A 2 


i/j ad (E,R)jp eS (R)dR 


(13) 


Since only l x — l y — L — S — 0 contributes, let us define 

Z K (r , R) = p _5/2 Xooooic(p)A(R-ph / / 2 2 ’ 1/2 (cos2«), (14) 

where Mr is a normalisation factor coming from Voooo/f gi yen by Eq. (10) of Ref. 0 and 
is a Jacobi polynomial 20]. After integration over all angles, the reduced transition 
probability m becomes 


POO 

V / u E (R)u§\R)dR 

K 


Rgx(E) — 6A^ 

It involves the Jl-dependent effective functions 

POO 

U^\R) = R Z K (r, R)u d (r)rdr , 

Jo 

the sum of which forms the radial part of ifj e g(R). 


(15) 


(16) 


III. RESULTS AND DISCUSSION 


A. Conditions of the calculation 


The central Minnesota interaction [2l| reproduces the deuteron binding energy and fairly 
approximates the low-energy nucleon-nucleon scattering properties. The deuteron wave 
function ip d is calculated over a Lagrange-Laguerre mesh involving 40 mesh points and a 
scaling parameter h = 0.25 fin (see Ref. jl3]] for details). An energy E d = —2.202 MeV is 
obtained. The calculations are done with h 2 /2mjv = 20.734 MeV fm 2 . 

The initial a + n + n bound state is calculated as explained in Ref. The number 
of components is limited to /i max = 24. The same nucleon-nucleon interaction is used, i.e., 
the Minnesota interaction with an exchange parameter u — 1. The a + N potential is 
however different from the one employed in Ref. 13]. Here we employ the a + N potential 
of Voronchev et al. m with a multiplicative factor 1.035 in order to reproduce the 6 He 
binding energy. This change of interaction is motivated by the fact that we want to use the 
same interaction for the derivation of the a + d folding potential. The renormalization factor 
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slightly affects the 5 He properties. For the 3/2 resonance, the original potential of Ref. ua 
provides an energy Er = 0.80 MeV and a width Y = 0.75 MeV, in nice agreement with 
experiment. Introducing the renormalization factor provides Er = 0.55 MeV and Y = 0.40 
MeV, but does not affect the unstable nature of the resonance. 

Since the valence neutron and proton in the 6 Li nucleus belong to the 0 p 3 / 2 subshcll, we 
use the p-wave a + N potential of Ref. flGj when deriving the a + d folding potential by 
using Eq. Q. For the s wave, this potential yields two bound states for 6 Li with energies 
E 0 = —19.87 MeV and E\ = —0.83 MeV, respectively. The first one is forbidden by the 
Pauli principle and the second one is underbound compared with the experimental ground- 
state energy E exp = —1.473 MeV. The a + N potential of Kanada et al. 0 employed in 
Ref. na does not yield an a + d folding potential with a physical bound state in the s wave. 

The numerical calculations of the Gamow-Teller matrix elements are done with several 
potentials, (i) The simple Gaussian attractive potential V a = —76.12 exp(—0.2r 2 ) 0 si¬ 
multaneously provides the correct 6 Li binding energy (together with a forbidden state) and 
a fair fit of the low-energy experimental phase shifts, (ii) The folding potential does not re¬ 
produce the 6 Li ground-state energy. Therefore, we multiply the central part of the original 
a + N potential by a factor fi = 1.068. The a + d folding potential Vf\ moves the physical 
state to E\ = —1.470 MeV. (iii) The folding potential Vf\ does not have the same quality 
of phase shift description as the simple Gaussian potential V a . Therefore, we also choose 
another multiplicative factor f -2 = 1.15 for the folding potential Vf 2 , which gives a stronger 
binding for the 6 Li ground state, E\ = —2.386 MeV. (iv) Finally, for the sake of comparison 
with Ref. jsi], we also perform a calculation with their Woods-Saxon repulsive potential (V))- 
Of course, it does not bind 6 Li. 

The s-wave phase shifts for the different a + d potentials are compared in Fig. Uwith the 
results of phase-shift analyses mm. The description of the a + d phase shift is poorest for 
the repulsive potential. Fair, almost identical, phase shifts are obtained with V a and V/ 2 - 

In Eqs. CHI) and (JHD, the integration over r is done by using the Gauss-Laguerre quadra¬ 
ture consistent with the Lagrange mesh. This ensures numerical convergence for the tran¬ 
sition probability. The integration over variable R in Eq. ca> is performed with the simple 
trapezoidal rule with a step 0.05 fm. Later we show that with this choice of step, convergent 
results for the transition probability are obtained with 600 points, which corresponds to a 
maximal a + d relative distance i? max = 30 fm. 

We have also calculated the Gamow-Teller matrix element for the f3 decay to the 6 Li 
ground state. For this calculation, we replace wave function in Eq. (fill) by the 

n + p wave function of the 6 Li ground state obtained with the same nuclear interactions 


a 


as for 6 He. The result Bq t = 4.489A 2 is about 5 % below the experimental value 7 ?gt = 
4.745A 2 . With the potential of Ref. [22] for the 6 He description, we find B GT = 4.636A 2 . 
The sensitivity with respect to the 6 He wave function is therefore small. 


B. Effective wave functions and their integrals 

In order to analyze the cancellation effects in the Gamow-Teller matrix element for the 
/3 delayed deuteron decay, we display in Fig. [3 the integrals 

Ie K \R)= [ u E (R')u { «\R’)dR' (17) 

Jo 

at the a + d relative energy E — 1 MeV for different K values. They are calculated by using 
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FIG. 1: s-wave phase shifts obtained with different a + d potentials: attractive Gaussian potential 
V a 0, folding potentials Vf\ and V /2 (see text), and repulsive Woods-Saxon potential V r 8]. 
Phase shifts are taken from analyses of experimental data in Refs. 0 (open dots) and |24| (full 
dots). 

the a + d potential of Ref. 0. The reduced transition probability is given by the limit 


■Bgt(E') — 6A“ 


lim J2i ( e\ R ) 

it —xxd ^—* 


K 


(18) 


^Frorri Fig. [21 one can see that at large R values the dominant contribution to Y2k ^e' f° r 
all K values up to K max comes from the K = 2 and K = 8 components. Components for 
K = 4 and K = 6 as well as for K > 12 are not visible with the linear scale of Fig. [21 
Although the K = 0 component is rather important around R — 5 fm, it is suppressed at 
large R values even more than the K = 10 component. 

To understand this interesting effect, we display different components u^\R) of the 
effective wave function as dotted lines in Fig. 01 The full line represents the sum 


R'u 


ua(R) = “£’(«) 


(19) 


K=0 


In Fig. 01 we also show the a+d scattering s-wave function ue for E — 1 MeV. It is important 
to note that this function keeps a constant sign in the interval 5-19 fm. This constant-sign 
interval is even broader for smaller values of E. The K = 0 and K = 2 components are 
dominant at all relative distances R. They exhibit a maximum below 5 fm. One observes 
that ufg keeps a constant sign over the whole region while u R changes sign at R ~ 2 fm. 

Since all other components are not visible in Fig. 01 we turn in Fig. 0to a logarithmic 
scale. For relative distances from R = 6 fm up to 25 fm, the contribution of the K = 8 
component is larger than the contributions of the K = 4 and 6 components. This is due to 
a zero around 10 fm in both components. 
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The curves in Figs. El and 0 indicate that the product ue(R)u^\r) for K = 0 changes 
sign several times. The integral I^\R) is first positive, starts decreasing at the first zero 
of ue, changes sign near 2 fin and increases again at the second zero of ue- The combined 
effect of both zeros results in a cancellation between the internal and external parts of the 
corresponding limit I^\+oo). These zeros at short distances are due to the existence of two 
bound states in the potential (one physical and one forbidden). The numbers and locations 
of zeros are typical of the 6 He /3 decay so that the cancellation should not occur for the 
decay of other halo nuclei. It would not occur so strongly with a single zero. 

For K = 2, the combined effects of the zero of u y e{ { and of the first zero of ue is just a 
small plateau near 2 fin in Fig. 0 The second zero of Ue gives a minimum near 5 fm. The 
resulting /Jy (+oo) for the K = 2 component also yields an important cancellation, but not 
as strong as in the K = 0 case. 

The effective functions for K = 4 and 6 are very small in the region where Ue keeps a 
constant sign and lead to negligible contributions. Let us recall here that the convergence 
of low -K components is not reached if Ji max is not large enough Q. 

A new situation appears for the K = 8 component. The effective wave function is much 
smaller than K = 0 or 2 but the cancellation is minimal since it does not change sign. 
Hence it gives the second largest I^ at infinity. The same mechanism applies for the 
smaller K > 10 components. The K = 10 integral still contributes significantly to the total 
sum. 

In Fig. 0 the integrals Ie(R) = ^e\^) calculated at the energy E = 1 MeV 

are represented for the different potentials. The repulsive potential V r displays a strongly 
different behavior from the other potentials. At this energy, V a and Vfi give almost the 
same result. This is due to the fact that, because of their similar phase shifts, their node 
near 5 fm in the scattering wave Ue is at nearly the same location. The comparison with 



R (fm) 

FIG. 2: Integrals I^\r) [Ea. H7|) ] at the energy E = 1 MeV for the a + d potential of Ref. 15] 
and different K values with A" max = 24 and i? max = 30 fm. The other components would hardly 
be visible at the scale of the figure. 
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FIG. 3: Effective wave functions u^\r) [Eq. (I16[) ] and u e g(R) [Eq. (1191) ] for the a + d potential of 
Ref. Q and different K values with K mSLX = 24 and if max = 30 fm. The scattering wave function 
ue at 1 MeV is also represented. 



FIG. 4: Same as in Fig. El in logarithmic scale. 

Vfi, where this node is about 1 fm farther away and leads therefore to weaker cancellation 
effects, shows the major role played by this node. 

The K = 2 and K = 8 components of the three-body hyperspherical wave function of 
the 6 He nucleus give dominant contributions to the integral Ie{R) at large values of R and 
thus to the Gamow-Teller reduced transition probability Bqj{E). This finding contradicts 
the assumption in Ref. [s|, that the K = 0 and 2 dominant contributions to the energy are 
sufficient to study this (3 decay mode. 
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FIG. 5: Integrals Ie(R) = YIk ^e^\R) at E = 1 MeV calculated with K max = 24 and R max = 30 
fm for different a + d potentials. 

C. Transition probability per time and energy units 
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FIG. 6: Transition probability per time and energy units dW/dE of the 6 He (3 decay into the a + d 
continuum calculated with the a + d potential V a and R max = 30 fm for several values of K max . 
Experimental data are from Refs. 5] (open dots) and |g] (full dots). The dotted line corresponds 
to the a + N potential of Ref. [ 22 I for the 6 He wave function. 

Before comparing transition probabilities with experiment, we discuss convergence as¬ 
pects. To study the convergence with respect to the value of the maximal hypermomentum 
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i f ma y, we display in Fig. 0 the transition probabilities [Eq. (fTTH) ] calculated with potential 
V a for if m ax = 2, 16, and 24. In each case the a + N interaction has been renormalized 
to reproduce the 6 He binding energy. In a logarithmic scale, the results for if max = 22 are 
essentially identical to those for if max = 24. From Fig. 0 we can see that using large if max 
values is crucial to obtain a good accuracy. This is analyzed in more detail below. 

Additionally, the convergence is faster for the repulsive potential V r and for the folding 
potentials Vj\ and Vf 2 - In the case of V r , the transition probabilities for if max = 16 and 
if max = 24 show almost the same features but have a larger value than for V a . However, 
even in this case, the choice if max = 2 in Ref. | 8 | is not realistic. 

A calculation with /f max = 24 performed with the original wave function of Ref. 0 , i.e. 
with the a + n potential of Kanada et al Hj, and the a + d potential V a is displayed as a 
dotted line. The results are somewhat larger and in less good agreement with experiment 
at low energies but they nicely reproduce both the shape and order of magnitude of the 
data. This indicates that the present model is not much sensitive to details of the model 
describing the 6 He wave function and confirms that convergence and an appropriate a + d 
potential are the crucial elements. 


TABLE I: Components I^\oo) as a function of K and If max at E = 1 MeV 


K A max = 2 

K 

Iv max 

= 16 

K 

il max 

= 18 

K 

Ji max 

= 20 

K 

Ji max 

= 22 

A max — 24 

0 3.67 x lO -2 

2.84 

X 

10~ 3 

4.02 

X 

10 -3 

5.21 

X 

10~ 3 

5.43 

X 

10~ 3 

5.74 

x 10 -3 

2 1.49 x 10” 1 

- 2.58 

X 

10~ 2 

- 2.09 

X 

10~ 2 

- 1.34 

X 

10~ 2 

- 1.17 

X 

10~ 2 

- 9.45 

x 10 -3 

4 

4.99 

X 

10~ 3 

5.27 

X 

10~ 3 

4.75 

X 

10~ 3 

4.90 

X 

10~ 3 

4.78 

x 10“ 3 

6 

1.14 

X 

10~ 3 

1.13 

X 

10~ 3 

1.56 

X 

10~ 4 

5.12 

X 

10~ 5 

- 3.06 

x 10“ 4 

8 

2.58 

X 

10~ 2 

2.66 

X 

10~ 2 

2.79 

X 

10~ 2 

2.84 

X 

10~ 2 

2.90 

x 10 -2 

10 

8.64 

X 

10~ 3 

8.90 

X 

10~ 3 

9.37 

X 

10~ 3 

9.55 

X 

10~ 3 

9.78 

x 10 -3 

12 

4.32 

X 

10~ 4 

4.15 

X 

10~ 4 

7.32 

X 

10~ 5 

8.38 

X 

10~ 6 

- 1.50 

x 10“ 4 

14 

2.03 

X 

10~ 3 

2.14 

X 

10~ 3 

2.21 

X 

10~ 3 

2.26 

X 

10~ 3 

2.30 

x 10 -3 

16 

2.01 

X 

10~ 3 

2.11 

X 

10~ 3 

2.33 

X 

10~ 3 

2.42 

X 

10~ 3 

2.53 

x 10 -3 

18 




4.82 

X 

10 -4 

4.70 

X 

10 -4 

4.71 

X 

10~ 4 

4.62 

x 10“ 4 

20 







1.49 

X 

10~ 4 

1.45 

X 

10 -4 

1.22 

x 10“ 4 

22 










3.26 

X 

10~ 4 

3.49 

x 10 -4 

24 













1.59 

x 10 -4 

sum 1.86 x 10” 1 

2.21 x 10" 2 

3.01 x 10“ 2 

3.92 x 10" 2 

4.23 x 10“ 2 

4.53 x 10“ 2 


The low value of dW/dE results from different cancellation effects which themselves are 
sensitive to the convergence of the different components of the three-body wave function. 
The final order of magnitude implies an accurate treatment of the convergence of the wave 
function and in particular of its halo part. In order to illustrate this mechanism, we display 
in Table [T| the numerical values of the components I^\oo) of the Gamow-Tcller integral for 
different values of the parameter K majX at 1 MeV. Let us start with /f max > 16. For each 
value of A' max , one observes that the dominant components are indeed K = 2 and K = 8 . 
Components beyond K = 16 become rather small. However, when K max is increased, 
all components I^\oo) of the matrix element are modified. As emphasized in Ref. [13. 
increasing A" max in the three-body model does not only mean adding components but, even 
more, improving the convergence of lower K components. This is illustrated when following 
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a row in Table HI The value of each component slowly converges when K max is increased. If 
the experimental data were much more accurate, higher values of K max should probably be 
considered. 

The dominant K = 2 and K = 8 components have opposite signs. This effect adds 
another level of cancellation in the matrix element. It increases the role of the other compo¬ 
nents and especially the collective role of high-A' components. Finally a comparison with the 
first column shows why a calculation restricted to K max = 2 has little meaning: the K = 2 
component is too large, has a wrong sign, and is not counterbalanced by other components. 
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FIG. 7: Transition probability per time and energy units dW/dE of the 6 He /3 decay into the 
a + d continuum calculated with the a + d potential V a and A' max = 24 for several values of A max . 
Experimental data are from Refs. _5| (open dots) and |fj (full dots). 


In Fig. [3 the transition probability dW/dE obtained with potential V a for a fixed A" max = 
24 is presented for different values of A max , i.e., 14 fin, 20 fm, and 30 fm. Calculations for 
Umax = 35 fm show that convergent results are obtained at A max ~ 30 fm. Taking properly 
account of the halo extension is very important in a correct treatment of the very small 
transition probability dW/dE of the 6 He (3 decay into a + d. 

In Fig. |HI we display the transition probability for different potentials, calculated with 
An iax = 24 and A max = 30 fm. The best description of the experimental data of Ref. [|j 
is obtained with the attractive potential V a . The worst results correspond to the repulsive 
potential V r) which has no bound state and for which the description of the s-wave phase 
shift at low energies is poor. With potential V r , the location of the nodes in the scattering 
wave does not lead to a cancellation (see Fig. ED- The folding potentials Vj\ and Vj 2 have 
intermediate behaviors. Potential Vf \ overestimates the recent data while potential Vy 2 pro¬ 
vides a better order of magnitude but its energy dependence disagrees with the experimental 
one. 

We have also performed a calculation with the original wave function of Ref. H3 and with 
a folding potential based on the a+n potential of Kanada et al 0 of Vf 2 type, i.e. multiplied 
by 1.30 in order to fairly reproduce the phase shift. The results are indistinguishable from 
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FIG. 8: Transition probability per time and energy units dW/dE of the 6 He /3 decay into the 
a + d continuum calculated with different a + d potentials for K rnaY = 24 and i? max = 30 fm. 
Experimental data are from Refs, fij (open dots) and @] (full dots). 


the curve labeled Vf 2 at the scale of the figure. The shape of this curve is thus mostly due 
to the node locations of the a + d wave function. 

The success of the deep Gaussian potential could be attributed to the fact that it simul¬ 
taneously reproduces both the 6 Li ground state binding energy and the s-wave phase shift 
at low energies. However the discussion of Figs. OH indicates that an important ingredient 
is the existence of two nodes in u e g. In order to test this assumption, we remove the non¬ 
physical ground state of V a by using a pair of supersymmetric transformations j2fj|. The 
resulting phase-equivalent potential V a 51 has exactly the same 6 Li ground-state energy and 
the same s-wave phase shift as V a but its scattering wave functions have one node less at 
small distances. The resulting dW/dE is about one order of magnitude larger and resembles 
the one obtained with the folding potential Vf\ (see Fig. 0). Notice however that Vfi has 
two bound states but does not well reproduce the phase shifts. A second phase-equivalent 
potential Vy" 2 is obtained by removing the b Li ground state from Vy' 1 with another pair of 
transformations. This repulsive potential has still exactly the same phase shifts as V a but 
no bound state. Its scattering wave functions have no node near the origin. As expected, 
the corresponding transition probability dW/dE is now very close to the one obtained with 
potential V r . The comparison emphasizes the crucial role played by the forbidden bound 
state, in addition to the physical 6 Li ground state, for reproducing the order of magnitude 
of the experimental data. 

The total transition probabilities for different potentials are given in the first row of Table 
UD The second row contains results corresponding to the experimental cutoff [&]. The values 
in the last columns are derived from experimental branching ratios and from the 6 He half 
life jb]. As expected from the previous discussion, the result obtained with the Gaussian 
potential V a falls within the experimental error bars of Ref. jjb]. The other results are too 
large, especially with the repulsive potential. 
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TABLE II: Total transition probability per second W (in 10 6 s x ) for the f3 decay of h He into 
a + d. 

_ Va Vfi Vf2 K Exp. [5] Exp. [6] 

E > 0 MeV 1.95 25.1 5.80 340 2.2 ± 1.1 

E > 0.37 MeV 1.54 18.2 3.32 246 7.6 ± 0.6 1.5 ± 0.8 


IV. A COMMENT ON 72-MATRIX FITS 


The 72-matrix method has been extended by Barker [11| to the (3 delayed deuteron emis¬ 
sion. It has been applied to analyze recent experimental results [§]. Like in other models, it 
is crucial in the 72-matrix method to take care of the l arg e extension of the halo. Without 
entering into details which are explained in Refs. (6j, II l| . 
external corrections proportional to the integral 


this is achieved by introducing 


Je(ci) = 


Uj(R) 

uAa) 


u B (R)dR , 


( 20 ) 


where a is the 72-matrix channel radius and ue is replaced by its asymptotic expression CD- 
The factor rq(a) eliminates the problem of normalizing the approximation for the initial 


11 


, the notation E % (72) is used for Ui(R)/ui(a). 


wave function Ui(R). In Ref. 

In the model of Ref. [rj3], the asymptotic form of the two-body a+dineutron system is 
employed for tq, 


K +2n (R) = exp [— ( 2 /i ( 1 2 )3 mjv | E B | / ft 2 ) 1/,2 72], ( 21 ) 

where Eb = —0.975 MeV. However, three-body asymptotics are rather different from this 
expression. In Eq. G5D> this role is played by the effective radial wave function n eff defined 
by Eq. (TBUl . In order to avoid the knowledge of three-body wave functions, we suggest here 
an expression, 


poo 

< +n+n (72) = 72 / p- 5/2 exp[-(27n N \E B \/h 2 ) 1/2 p\ru d (r)dr, (22) 

Jo 

which is the projection of three-body asymptotics jTol. Ell ] on the deuteron wave function. 
For a pointlike deuteron described with Ud(r) oc r~ 1 S(r ), this function becomes 

< +2n ’ cor-(72) = 72~ 3//2 exp [— (2/r (12)3 mjv | E B \ / ftjy^R]. (23) 

It differs from Eq. m by the power factor 72 3//2 . This simple expression also deserves 
being evaluated. 

In Table mU we present the integral (HD calculated with < +2n , < +n+n , M “+ 2 «- cor - ; and 
u e g at two typical energies for different values of the channel radius a. One observes that 
the results obtained with the two-body asymptotic expression m are rather far from the 
realistic values obtained with u e g, even for a = 6 fin. A much better approximation is given 
by the three-body asymptotic expression (1221) . especially at higher relative energies. The 
corrected two-body approximation is smaller than the three-body approximation and not 
really close to the reference results. 
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TABLE III: External integrals [Eq. (120H ] in ii-matrix fits: X^ +2n , x ^ 2n ' cor ', 2^ +n+n , and 2 J? 
are calculated with uf +2n [Eq. iffll) ]. u “+ 2n ’ cor - [Eq. 031) ]. u^ +n+n [Eq. 021) ]. and u e g [Eq. (fTTH) ] . 
respectively, for different values of the channel radius a (in fm) and the relative energy E (in MeV). 


E a 

q-a.+2n 

/T-a+2n,cor. 

'■7-a+n+n 

X E 

7"eff 

X E 

0.5 4.0 

1.162 

0.224 

0.470 

0.598 

4.5 

1.391 

0.384 

0.651 

0.751 

5.0 

1.616 

0.552 

0.834 

0.913 

5.5 

1.834 

0.726 

1.018 

1.082 

6.0 

2.046 

0.903 

1.201 

1.432 

1.0 4.0 

1.464 

0.367 

0.698 

0.882 

4.5 

1.754 

0.601 

0.952 

1.091 

5.0 

2.024 

0.839 

1.200 

1.307 

5.5 

2.172 

1.075 

1.437 

1.522 

6.0 

2.496 

1.303 

1.660 

1.732 


In i?-matrix theory however, because parameters are fitted, the absolute normalization of 
the integrals displayed in Table ITTT1 is not important. It can be absorbed in a renormalization 
of the constants. The crucial property is their energy dependence. The four types of external 
integrals are shown in Fig. |H1 as a function of energy for the typical value a = 5 fm. One 
observes that x“+ 2n > cor - ; 2^ +n+n , and Zjf have very similar energy dependences. On the 
contrary, Z^ +2n displays a different shape. Hence, using this approximation in /?-matrix fits 
may signihcantly distort the energy shape of the f3 delayed deuteron spectrum. The a+7i + n 
approximation offers a very good approximation of the model results. Nevertheless, the 
corrected two-body expression provides the simplest significant improvement for i?-matrix 
calculations. 



FIG. 9: External integrals X^ +2n , x“ +2n,cor ', 2^ +n+n , and 2 J? as a function of the relative energy 
E for a = 5 fm. 
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V. CONCLUSIONS 


In the present work, we studied the /9-decay process of the 6 He halo nucleus into the a + d 
continuum in the framework of a three-body model. Three-body hyperspherical bound-state 
wave functions on a Lagrange mesh and two-body a + d scattering wave functions have 
been used. For the calculation of the /9-decay transition probabilities per time and energy 
units, several a + d potentials were tested: an attractive Gaussian potential 0 with a 
deep forbidden bound state, folding potentials derived from the a + N p -wave potential of 
Ref. [16], and a repulsive potential jjj|. 


We confirm that the low experimental values result from a strong cancellation in the 
Gamow-Teller matrix element describing the transition to the continuum ]9], This cancel¬ 
lation occurs between the internal and halo parts of the matrix element and is thus very 
sensitive to the halo description. Reaching convergence is not easy: the two-body and three- 
body wave functions must extend up to about 30 fm. From the analysis of the theoretical 
results, we have found that converged results require large values of the maximal hyper¬ 
momentum iF max . The dominant contributions to the transition probability come from the 
K = 2, K = 8, and K = 10 components of the three-body hyperspherical wave function. 
The K = 0 contribution is small due to an almost perfect cancellation of its internal and 
external parts in the Gamow-Teller matrix element. The K = 2 and K = 8 components have 
opposite signs which enhances the importance of other, and especially high-JF, components. 

The experimental transition probabilities per time and energy units of Ref. (!] are well 
described with the deep Gaussian potential of Ref. 0 which fairly reproduces the 6 Li 
binding energy and the a + d s -wave phase shifts. The quality of the agreement arises from 
the node structure of the initial and final wave functions in the internal part. With the 
help of phase-equivalent potentials derived with supersymmetric transformations, we have 
shown that the role of the forbidden state is also essential. We realize that the efficiency 
of the deep potential may be somewhat fortuitous since the nodes of the scattering wave 
function have to be at very precise locations. The fact that the data can be reproduced does 
not mean that the present model or the simple Gaussian potential are perfect. However 
the reasonable agreement with the data obtained with the same a + d potential but another 
6 He wave function indicates that the present model interpretation should be trustable. Most 
importantly the existence of a good agreement with experiment points toward the ingredients 
that are crucial in the interpretation of the /3 delayed deuteron decay of 6 He. One can expect 
a completely different behavior for this /9-decay mode in the case of 11 Li. Indeed the 6 He-case 
cancellations require precise numbers of nodes and precise locations of these nodes and it is 
very unlikely that this could occur so perfectly in another case. 

Our results allow testing the validity of external corrections necessary in the /?-matrix 
method jll]. We have shown that, in order to reduce a systematic bias in the integrals 
over the external region, the two-body asymptotics can usefully be replaced by three-body 
asymptotics or, more simply, be corrected by a factor f?~ 3 / 2 . 

Further progress on this /9-decay mode must come from consistent fully microscopic de¬ 
scriptions of the bound and scattering states. The results obtained with a microscopic cluster 
model still agree qualitatively with the most recent data [1] but overestimate them by 
about a factor of two. Progress may be expected from the possibility of calculating 6 He wave 
functions ab initio f§7j from realistic two- and three-body forces. However the present study 
shows that a successfull description of the /3 delayed deuteron emission will require very 
accurate bound-state wave functions up to distances as large as 30 fm and a development of 
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consistent scattering wave functions. The accidental cancellation occurring in this process 
will make a fully ab initio description particularly difficult. 
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